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Abstract 

A Deformed Gaussian Orthogonal Ensemble (DGOE) which interpolates be- 
tween the Gaussian Orthogonal Ensemble and a Poissonian Ensemble is con- 
structed. This new ensemble is then applied to the analysis of the chaotic prop- 
erties of the low lying collective states of nuclei described by the Interacting 
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Boson Model (IBM). This model undergoes a transition order-chaos-order from 
the SU{3) limit to the 0(6) limit. Our analysis shows that the quantum fluctua- 
tions of the IBM Hamiltonian, both of the spectrum and the eigenvectors, follow 
the expected behaviour predicted by the DGOE when one goes from one limit 
to the other. 

It is widely assumed that Random Matrix Theories (RMT) provide a basis to 
study quantum chaotic systems. In particular, it is expected that fluctuation properties 
of fully chaotic systems with time reversal symmetry follow the Gaussian Orthogonal 
Ensemble (GOE) whereas non-chaotic ones follow the Poissonian Ensemble. Some 
physical systems however may exhibit statistics intermediate between these two limits, 
as recent investigations have shown in the case of the excitation spectra and intensities 
of deformed 00 and spherical [^] nuclei. The analysis that has been performed on 
these systems have a more or less empirical character in the sense that they are not 
based directly on an ensemble of RMT. It is therefore important to test the reliability 
of the RMT predictions in these intermediate situations. 

Recently, two of us have constructed a Deformed Gaussian Orthogonal Ensemble 
(DGOE) using the maximum entropy principle applied to generic random matrices 
subjected to appropriate constraints^. The constraints imposed were such that the 
ensemble obtained goes from a pure GOE to a combination of two GOE's. This corre- 
sponds to the case of SU (2) symmetry breaking, when a quantum number which can 
have only two values, e.g. isospin, is partially conserved^]. To deal with the more 
general case of the GOE-Poisson transition an extension of the DGOE of Ref.|^ is re- 
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quired. In this paper we provide this extended DGOE and apply it to the analysis of the 
transition chaos-order exhibited by the Interacting Boson Model (IBM) Hamiltonian. 

We start using projection operators Pi =\ i >< i \, Qi = 1 — Pi, where | i >, with 
i = 1,2, ■..,N, are N abstract basis vectors, to split a generic Hamiltonian operator H 
in its diagonal and off-diagonal parts, Hq and Hi. Namely, 



N 
1=1 



N 

Hi = J2P,HQ, (2) 

i=l 

where Hq and Hi satisfy the identity H = Hq + Hi. 

Following the discussion of Ref. we add to the usual GOE constraints 

<TrH^>= J dHP(H)TrH^ = 1^ (3) 
< 1 >= y dHP{H) = 1 (4) 

the additional one 



TV 

<TrH^,>=J2 / dHP{H)Tr{PiHQiHP,) = iy. (5) 

i=i 

Maximizing then the entropy subjected to the above conditions we get, after nor- 
malization, the probability distribution: 



3 



P{H) = PGOE{H)exp[~PY.TrimQ,HP,)]{l + (^"^ (6) 

i « 

where a and j3 are Lagrange multipliers and 

Pgoe{H) = 2-^/2 (^— j exp(-a Tr Z/^) (7) 

is the GOE distribution. 

It is clear from Eq. |^ that when /3 — >^ we recover the GOE distribution while with 
the increase of (3 the system becomes less chaotic; in the limit /3 ^ 00 we recover the 
Poissonian Ensemble . 

Using an idea suggested by Dyson [0, Alhassid and Levine0 discussed an interme- 
diate ensemble as a solution of a non-equilibrium problem defined by an appropriate 
stochastic Langevin equation. The non-integrability is characterized by a parameter e 
which is the ratio between the variances of the non-diagonal to the diagonal elements 
of H. We easily find the relation e = , ^ showing that e goes from zero to unity as 
j3 varies from infinity to zero. 

At this point, we remark that by using the projectors Pi and Qi we can alternatively 
decompose the Hamiltonian H as 

N 

H = Ho + eJ2 PiHcQr (8) 

i=l 

where He = H{(3 = 0) showing that the problem may be reformulated in such a way 
that the parameter e appears as a coupling constant. This relation makes a connection 
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of our ensemble to the one discussed recently by Lenz and Haake0. 

We present now the numerical results. In Ref. a new parameter e was introduced 
such that the GOE limit is expected to be approached when e 1. This parameter 
may be expressed as 



Fig.l the parameter uj of the Brody distribution that fits the level spacing as a function 
of the size N of the matrix, keeping e or e fixed. We fixed the value of e to be that 
at N=50, the lowest dimension in the calculation. It is clear that lo saturates quite 
rapidly when e is kept constant. 

In Fig. 2, we have considered matrices of size N=200 and calculated the level spacing 
P(s), the spectral rigidity A3(L) and the distribution P(ln?/), where y is the square of 
the component normalized with respect to its average {y = (? j < >), for several 
values of the parameter e. We see that indeed we are very close to GOE when e = 1. 
The level spacing distribution exhibits the universality law which says that the level 
repulsion only disappears in the Poisson limit, e = or /3 = oo. We remark also that 
our P(s) are practically identical to those obtained by Lenz and Haake with N = 500 
matrices. With respect to P(ln?/), the distributions get broader as e decreases but we 
cannot fit them with only one Py. Actually, as in the previous case of the DGOE for 
two GOE's, an excellent fit is obtained if we use instead two P^'s. This behavior may 
be understood as a consequence of the fact that in the transition from the GOE to the 




(9) 



and for N=2 e = e while for large N, e ~ \fNt. To test this idea we have plotted in 
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Poissonian case, the components of the eigenvectors go from a situation in which they 
are equally distributed among the basis states, to a limit situation in which only one 
intensity, say Ui, is different from zero. Near this limit, the average of the components 
becomes equal to < y >~ 1 — ?/j and this explains why the distribution shift to the 
left in Fig. 2 . This fact, together with the splitting of the intensities in two sets fitted 
by two different Pjy's, seem to be universal features of the eigenvector statistics in the 
transition from a chaotic to a regular regime. 

The Interacting Boson Model (IBM) has been successfully applied to the phe- 



nomenological description of low-lying collective states of atomic nuclei |T^. In Ref. [g] 
the fluctuation properties of such states were analyzed in the framework of IBM. We 
turn now to a detailed analysis of the statistical properties of the IBM Hamiltonian, 
using the above intermediate ensemble. The IBM Hamiltonian may be written as 

H = rind-{l-7])Q^-Q'' (10) 
where is the number of d bosons and is the quadrupole operator 
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S-d (111 



Qx = (^t X s + sT X J)(2) + x{d^ X J)(2), (12) 

with the relation dfj_ = {—^d-^. The six bosons s and d^ span a six- dimensional 
Hilbert space which has U{Q) symmetry. The IBM Hamiltonian has three dynamical 



symmetry limits which describe vibrational nuclei for rj = 1, rotational nuclei for 
7] = and X = ~a/7/2, and, finally, 7-unstable nuclei for i] = and x = 0- It was 
shown recently that among these three limiting cases, there is a region where the 
Hamiltonian becomes chaotic. In particular, taking rj = and various values of x in 
the range —\/7 /2 < x < 0, namely going from rotational nuclei to 7-unstable ones we 
encounter the following results shown in Fig.3. With regards to the level spacing and 
the A3 we are just reproducing the results of Ref. [|]. However, instead of presenting 
the BE(2) distributions as was done in Ref.p|, we show in the first column of the figure 
the statistics of the components. This is necessary in order to compare with the DGOE 
statistics presented in Fig.2. Besides, the components were calculated with respect to 
the basis of the nearest regular cases. 

The IBM spacings and A3 distributions show a clear Poisson-GOE-Poisson transi- 
tion as one varies x from -1.25 to -0.1. The two transitions considered separately have 
a behavior which compares well with the description given by the ensemble. However, 
insofar as the component distributions are concerned, there is a marked difference, 
which seems to arise from what we may call the intermediate structure that modulates 
the components distribution. This intermediate structure was refered to as secular 
variations in Ref.p|. Therefore in order to make a sensible comparison we still have 
to extract from the components this secular variation. This can be performed by in- 
troducing a local average defined by an appropriate moving window. Following the 
procedure of Ref. we have used a Gaussian window and have taken the labels of the 
components as variables. Explicitly, we have used 
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and the normalized components are defined as = • In fig-4 we show the 

resulting distributions. In the column on the left the we have the distributions for 
the DGOE in the case of 200x200 matrices whereas the other two present the IBM 
results. The remaining two columns show the transitions order-chaos starting near the 
rotational limit {x = —1.25, bottom) and starting near the 7— unstable limit (x = 
—0.10, bottom), respectively. We can see that the distributions for the two transitions 
of the IBM Hamiltonian follow the pattern described by the ensemble. 

In conclusion we have shown in this paper that the extended version of the DGOE 
of Ref. 1^, can indeed describe well the statistical behavior of a realistic nuclear model 
such as the IBM. Similar observations can be made with regard to the study of high- 
spin states 0, and the analysis of spherical nuclei made in Ref.0]. We have reasons 
to believe that any physical system is said to exhibit order-chaos transition when its 
statistics follows rigorously that of a extended DGOE. 
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Figure Captions: 

Fig.l-The Brody parameter, cu, of the adjusted spacing distributions, vs. tlie size 
of the matrix. The diamonds correspond to e = 0.01, whereas the crosses to e = 
e/VF^ = 0.07/ ViV^. 

Fig.2-The distributions, P(ln?/), P(s), and A3 of the DGOE for N=200 and various 
values of e. See text for details. 

Fig.3-Same as Fig. 2, for J=6+ IBM states. See text and Ref. for details. 

Fig.4-The components distributions after normalization, the Ihs column presents 
the DGOE 's and the other two the IBM ones . See text for details. 
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